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Abstract 



We study boundary bound states using the Bethe ansatz formalism for the 
finite XXZ (A > 1) chain in a boundary magnetic field h. Boundary bound 
states are represented by the "boundary strings" similar to those described 
in Q. We find that for certain values of h the ground state wave function 
contains boundary strings, and from this infer the existence of two "critical" 
fields in agreement with Q. An expression for the vacuum surface energy 
in the thermodynamic limit is derived and found to be an analytic function 
of h. We argue that boundary excitations appear only in pairs with "bulk" 
excitations or with boundary excitations at the other end of the chain. We 
mainly discuss the case where the magnetic fields at the left and the right 
boundaries are antiparallel, but we also comment on the case of the parallel 
fields. In the Ising (A = oo) and isotropic (A = 1) limits our results agree 
with those previously known. 
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One-dimensional integrable quantum field theories with boundary interactions |J have 
been intensively studied recently because of their applications in condensed matter physics 
(see e.g., Q). A powerful method for dealing with such problems is the Bethe ansatz, 
allowing basic physical properties to be extracted from the system of coupled transcendental 
equations. Among others, it allows to solve the boundary sine-Gordon model via its lattice 
regularization, the inhomogeneous XX Z (|A| < 1) chain in a boundary magnetic field |E], 

In this Letter we study the XXZ chain with even number of spins L in a boundary 
magnetic field, 

H = \ + + + h*>i + ^ J . (!) 

in the regime A > I, h± > 0, h 2 < 0, focusing on the effects peculiar to systems with 
boundaries 0. At h\ — hi — this model describes one-dimensional antiferromagnet with 
non-magnetic impurities, accessible experimentally. We exploit the Bethe ansatz solution for 
this model, first derived in |7[], together with the well-known results for the periodic chain ||. 
We find new "boundary string" solutions to the Bethe equations, similar to the boundary 
strings existing in the |A| < 1 regime For certain values of the boundary magnetic 
field the ground state configuration contains boundary 1-strings. Boundary excitations are 
obtained by removing (or adding, depending on the sign of h) boundary strings from the 
ground state wave-function. Their energy was first obtained in by the algebraic approach. 

A peculiar feature of the Bethe ansatz solution of the periodic chain is that the excitations 
(holes in the Dirac sea) appear only in pairs ||. We argue that similarly the boundary 
excitations can appear only in pairs with bulk excitations or with boundary excitations at 
the other end of the spin chain. This is not evident in the algebraic approach of |2j . 

Using the Bethe ansatz solution we calculate the surface energy (see e.g., [fLOfl): 

E sur f(L, A, h) = E gr — Eg r , (2) 



in the thermodynamic limit L = oo. Here E gr is the ground state energy of ([!]) and Eg r is 
that of the periodic chain. We give an interpretation of our results in the limits A — > oo and 
A — > 1, corresponding to the ID Ising and XXX models respectively. Finally, we comment 
on the structure of the ground state when the boundary magnetic fields are parallel. 

Let us first set up the Bethe ansatz (BA) notations and list the relevant results about 
the XXZ chain 0J|. In the eigenstates of (0) were constructed for arbitrary A. As 
usual in the BA picture, the n-magnon eigenstates \n), satisfying 7i\n) = E\n), are linear 
combinations of the states with n spins down, located at sites x\, ...,x n : 



contains n parameters pj G (0, 7r) which are subject to quantization conditions, called Bethe 
equations (BE): 




The wave-function 




(3) 



p 



e ip J +hx-A e ip i + h 2 - A 



n 



e 2iL Pj _ 



1 + (hi - A)ew 1 + (h 2 - A)e^ 



I J e **(Pi,Pi)_ 



(4) 



The energy and spin of the n-magnon state are given by 0: 




(5) 



It is convenient to rewrite BE using the following mappings: 



A = cosh 7 > 1, 



7>0, 



(6) 




(our definition of p(a) differs from that of by the shift a — > a + tc and it was chosen in 
such a way that p(a) be an odd function that maps — tt < a < n to — 7r < p < 7r), 
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sinhm-H) 7/ 
/i = cosh7 + — —— — — — = sinh 7 ■ coth -(if + 1), h Hm < \h\ < oo, (8) 

sinh ^(1 + if) 2 

cosh |(1 -If) i 

/i = cosh7 = smh 7 • tanh - (if + 1), \h \< h Um . (9) 

cosh ^(1 + if) 2 

The latter two mappings are defined on if G (—00, 00) and are necessary to cover the region 
—00 < h < 00, with positive h corresponding to if G (— l,oo). The value hu m = h(oo) = 
sinh 7 lies between two critical fields ha/ , h& defined as follows @ : 

hg}=A-l, h®=A + l. (10) 

Both critical fields correspond to if = 0, and the gap < h < hhfi corresponds to 
< H < 00. In these notations eq. (f|) becomes: 

cosh \ {i(Xj + 7) ] 2L ^ if ) TT ~ ^ m S ^ n ^ ^ m 



cosh \{%a.j — 7) 



sinh ^iay — za m — 27) sinh |(iaj + ia m — 27) 

(11) 



where 



cosh 1 (za + 7 if ) , , 

# = u 1 • 4k> ^™ < /1 < oo, (12) 

cosh ^(za — 7-n ) 

sinh i (ia + jH ) , , 

S a, If) = — A < h lim , (13) 

smh ^(za — 7-n ) 

are called boundary terms. The energy eq. @ takes the form: 

1 n inh 

E = -{(L - 1) cosh 1 + h 1 + h 2 ] - 2 sinh 7 V j/(a) = — . (14) 

2 ^-^ cosh 7 + cos a 

In the thermodynamic limit L — > 00 the real roots a, of BE form a dense distribution in 
the open interval (0,7r) with density p(a), c£f = 2L(p + p h )da being the number of roots in 
the interval da. The logarithm of eq. ( |iT| ) is: 
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11 n 

2Lp(ctj) + - In B(atj, Hi) + - In B(aj, H 2 ) + 4>{2aj) = (f>(aj — a/) + 4>{ a j + a i) + 27r/j, 

1=1 

(15) 



where form an increasing sequence of positive integers, and 

^ = - iln ^Wi^y *<°) = a (16 > 

Taking the derivative of eq. ([□]) and defining p for negative a by p(a) = p(—a), we obtain 



P'(a) + ^rp' bdry (a)= I (f>'(a- P)p(P)dp + 27r(p(a)+p h (a)), (17) 



with 

y B[a,Hi) B(a,H 2 ) 

The presence of delta-functions in fll8|) is due to the fact that ctj = and a, = n are always 
solutions to flllD , which should be excluded, since they make the wave- function (|3]) vanish 
identically 0. 

In eq. (|T7|) the "boundary terms" are down by the factor 1/2L. Neglecting p' bdry and 
setting ph = 0, we obtain the equation for the ground state density of the periodic XXZ 
chain JTTJ] . Solving it by the Fourier expansion 

/(«) = E /(0e a °, /(0 = ^ T f(a)e- Ua da, (19) 
and using flI3|), we recover the result for the ground state energy of the periodic chain jll 



2np per (n)= P[H - , 4>\n) = e- 2 ^, p\n) = (-l)V^H (20) 
1 + 0'(n) 

LA f 11 LA 00 ^ e ~7l«l 

^gr = ~s 2Lsinh7 / p per (a)p'(a)da = — Lsinh7 — . (21) 

y 2 /_ 2 ' cosh 771 

n=— oo 

The spin of the ground state is S z = L/2 — L f_ p per da = 0, which is the well-known result 

hi. 
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An elementary "bulk" excitation above the vacuum in the model ([1]) is a hole in the 
distribution of Ij, but only a pair of holes can occur for the periodic chain, as argued in 
0. Thus physical excitations contain an even number of holes. The energy of the hole with 
rapidity 6 can be easily computed: 

e h {d) = sinh 7 V ^4 > 0, (22) 

^— ' cosh 7n 



and the spin with respect to the vacuum is S z = 1/2. (Our result, eq. fl22"|) , differs from the 
conventional one by the shift 9 — > 6 + n, but the dispersion relation is unchanged by rapidity 
reparametrization.) 

Analogous arguments can be applied to analyze "bulk" string solutions with complex 
a. Although there exists an infinite hierarchy of complex strings of arbitrary length and 
quartets, their energy vanishes with respect to the vacuum fl2| . 



So far we discussed the bulk excitations, which are essentially the same as in the periodic 
chain. Let us now turn to the new solutions of eq. (0), boundary strings. The analysis is 
close to that of Boundary excitations have their wave-function (§) localized at x = or 
x = L, and in the limit L — > 00 the two boundaries may be considered separately. Let us 
study first the left boundary, h\ > 0. The fundamental boundary 1-string consists of one 
root located at ao = —ijHi for < h\ < , and at «o = tt — i^Hx for har < h\ < 00 
(in both cases — 1 < H\ < 0). It is a solution of BE due to the mutual cancellation of 
the decreasing modulus of the first term in (|TT| ) and the increasing modulus of the second 
term B(a, Hi) as L —>■ 00 and a — > oq. When < h x < hf) , no such solution exists. 
Introduction of such a string into the vacuum whose root density p(a) satisfies 

p'{oc) + ^ P ' bdry (a) = J 4>'{a - (3)p((3)df3 + 2vrp(«), (23) 

leads to the redistribution of roots described by 5p = 2L(p — p). 5p is a solution of the 
equation 
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= / <f)'(a- (3)5p(f3)d(3 + (f)'(a-a ) + (j) > (a + a ) + 27T6p, (24) 

J — 7T 

Solving (|24D , we find 

2cosna e- 2 ^l 

27r,5p(n) = ^ - ^ H . (25) 

For the energy and spin of the boundary 1-string with respect to the Dirac sea we obtain: 

5pp'(a)da = — sinh7 > , S z = —1/2. 

_ ' cosh7n 

n=— oo 

(26) 

We see that the energy ( ^6|) is negative, so the correct ground state wave-function (|3]) should 
contain the boundary 1-string root ao when the value of hi is not in the gap h& < hi < h& . 
The ground state density p in this case satisfies the equation 

p'(a) + ^(jp'uryia) - <P'(a - a ) - <P'(a + a )) = J <P'(a - f3)p{f3)d{3 + 2^p{a). (27) 

The boundary excitation is obtained by removing from vacuum the root cto, which means 
that it has the energy —it > and spin 1/2, equal to the spin of the bulk elementary 
excitation. Substituting the value of a into (ffiSl), we get the boundary excitation energy, 
which precisely agrees with the one obtained in 0: 

e b {hi) = sinh 7 V ^-L- , -KHi< 0, (28) 

*— ' cosh 771 

n=— oo 

with k = 1 if hi < h^cr and k = 2 if hi > h% . 

From ( p2| ) and (p8|) we see that for hi < h$ the energy of the boundary excitation is 
smaller than the energy of the zero-rapidity bulk excitation, and becomes equal to it at 
hi = her ■ So in this regime we may interpret the boundary excitation as the bound state of 
the kink, which gets unbound at hi = h^r ■ For hi > h!$ such interpretation is impossible 
because the energy gap is smaller than the boundary excitation energy. 



Besides the fundamental boundary 1-string, there exists an infinite set of "long" boundary 
strings, consisting of roots o,q — 2i&7, «o — 2i(k — 1)7, «o + 2m7 with n, k > 0. We will call 
such solution an (n,k) boundary string (thus the fundamental boundary string considered 
above is the (0,0) string). One can use the same arguments as given in Q to show that the 
(n,k) string is a solution of BE when its "center of mass" has positive imaginary part and 
the lowest root «o — 2z&7 lies below the real axis. Thus sufficiently long boundary string 
solutions exist even in the region h$ < hi < hcP . However, a direct calculation shows 
that their energy vanishes with respect to the vacuum, so they represent charged vacua. [] 



(Analogous phenomenon occurs for the "long" strings in the bulk []12~ |: if the imaginary part 
of a lies outside the strip —27 < Ima < 27, it gives no contribution to the energy.) For 
< hi < h^cr and h& < hi < 00 the (n,0) strings also represent charged vacua, while (n,k) 
strings with k > 1 have the same energy (ESp as the boundary bound state found above, and 
hence represent charged boundary excitations. [] 

Consider now the right boundary, hi < {H 2 < —1). Now the fundamental boundary 
1-string solution = —i r )H 2 exists for any value of h 2 in the interval —hi im < h 2 < (resp. 
Oq = tt — for h 2 < —hu m ). Explicit calculation shows that it has non- vanishing energy 
only if —2 < H 2 < —1, which corresponds to —h$ < h 2 < (resp. h 2 < —h^r)- For such 
values of h 2 the energy of the 1-string with respect to the Dirac sea is positive and equal to 
£b(—h 2 ) (see eq. (p8|)), and its spin is S z = —1/2. In some sense the pictures are dual for 
the positive and negative h cases: there exist two states when \h\ is not between \h^)\ and 



1 Ks an example, consider the boundary (l,0)-string consisting of the roots ao + 2^7, ocq. It exists 
if — 1 < Hi < 1, although the (0,0)-string exists only if — 1 < Hi < 0. The (l,0)-string has charge 
S z = — 1 and vanishing energy with respect to the vacuum. 

2 E.g., (1,1) string with roots ao + 2z7, oo, «o — 2ry has S z = —3/2 and energy given by (28) with 
respect to the physical vacuum. 
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(2) 

\hcr \i one with a boundary 1-string and one without. One of them is the ground state and 
another is the excited state at the boundary, and these states exchange their roles when the 
sign of h changes. The analysis of long boundary strings is very similar to that at the left 
boundary, and therefore will be omitted. The net result is again that long boundary strings 
represent charged vacua or charged boundary excitations. 

In all examples shown above, the charge of boundary excitations turned out to be half- 
integer. One can easily check that this is true for all boundary strings representing charged 
excitations. Since the charge of physical excitations is obviously restricted to be an integer 
(see (H)), we conclude that a boundary excitation can appear only paired with the bulk 
excitation of half- integer charge (i.e. containing an odd number of holes), or with a boundary 
excitation at the other end of the chain. We give a qualitative interpretation of this fact 
below. 

To compute the vacuum surface energy, eq. (Q), of the model ([]]), one should use eq. 
([14]) in the limit L = oo with the root density determined from eqs. (|23| ) or ( p?D and the 



boundary terms (|12|) or (|13D , depending on the value of h. Define for convenience 



^ = - + 23^^^-. (29) 

n=l 

We consider separately the following intervals for positive hi and negative h 2 : 

1) |^i,2| < her . The ground state contains one boundary 1-string, corresponding to h\. 
The spin of the ground state can be found to be S z = 0. Using eqs. (|T3|) , (|18|) , and (|27|) and 
subtracting the bulk contribution ( pl| ) we get 



E surf = -(h + h 2 ) - g(A) - sinh 7 V(-l) n - ^ . (30) 

2 ^-^ cosli7n 

71=1 ' 

(2) 

2) |^i,2| > her . The ground state contains one boundary 1-string and has S z = 0. From 
eqs. flT2|) and (0) it follows: 

E surf = A h i + h) - p(A) - sinh 7 V ^ . (31) 

2 ' cosh7n 

n=l 
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3) h cr < < hum- The ground state has no boundary strings and its spin is zero. 



From ( |23| ) and ( p~3|) one obtains the same expression as in case 1). 

4) hu m < \hi )2 \ < lihr ■ From fl2"3] ) and fli"2] ) one obtains the same expression as in case 2). 



The ground state has the same structure as in case 3). 

A qualitative plot of the surface energy as a function of h {h = h\ — — /12) is given in 

Fig.l. The apparent difference between (^) and ([H]) is an artefact of our parametrization 

of h in terms of H. In fact, E sur f is an analytic function of h in the domain h £ (0, 00), 

which can be seen after substituting if as a function of h according to In this sense 

(1 2) 

the fields hcr are not actually "critical." We find for hi = h 2 = the value 

2-^1 \ _|_ e 4n7 " 1 ~ 2-^ 
n=l n=l 

In the extreme anisotropic limit A — > 00, h ~ A of the XXZ chain (JJ) one gets the 



E 



surf 



A / 1 ^e 2 " 7 -l ^ _i n \ 

— + 4snh 7 -+> — 5— . 32 



one-dimensional Ising model: 

'L-l 



w = J 1 E Aa « <^ + ^ + ^ > ( 33 ) 



i=l 



2 

In this limit from (§)-(§) one has 

h « A ± e~ 7 ", (34) 

and the gap between /li* and fo£r dissapears, so for any h there exists a boundary bound 
state. The energy of the "bulk" hole ( p2|) becomes independent of 6 and equal to A, since 
only n = term contributes to the sum. The energy of the boundary bound state ( f28|) 
becomes £& = A ± e _Tf/ = /i. This suggests the following interpretation in terms of the Ising 
chain: in the Ising ground state the i-th spin has the value (— 1)\ Local bulk excitation 
of the smallest energy 2 A can be obtained by flipping one spin (the first spin excepted). 
The arising two surfaces (domain walls) separating the flipped spin from its right and left 
neighbours are called kinks and carry the energy A each. The kink corresponds to a hole 
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in the Dirac sea in the Bethe ansatz picture, and these holes can exist only in pairs, which 
is obviously the case for the kinks in the Ising chain. The charge of the one-spin-flipped 
state is equal to one, in agreement with the charge of two holes in BA. In addition to charge 
1 excitation, one has charge excitation of the same energy obtained by flipping any even 
number of spins in a row. In the BA this corresponds to the "2 holes and 2-string" state. The 
boundary bound state is obtained by flipping the first (or the last) spin. Such a state has the 
energy h% + A above the vacuum energy in the Ising limit, where h\ is the contribution of the 



boundary term in eq. fl33|) and A is the energy of the kink created because of the interaction 
of the boundary spin with the bulk ones. Thus flipping the boundary spin actually gives a 
combination of the boundary excitation and the bulk kink. Still another possibility is to flip 
all spins, creating two boundary bound states, one at each boundary. This explains why, in 
the BA picture, a boundary excitation can exist only if paired with a hole in the Dirac sea 
or with another boundary excitation. 

The ground state surface energy (§) of the Ising chain in the thermodynamic limit is 
(A — h\ + h 2 )/2. The A/2 contribution here is the bulk interaction energy that we lost 
when we disconnected the periodic chain, and ±/i 12 /2 is the contribution of each of the 
boundary terms. Taking the limit 7 —>■ 00 in eqs. (j30|)- (j3l|), we obtain the expected result 
(A - h x + h 2 )/2. 

In the isotropic (rational) limit A — > 1 (i.e., 7 — > 0) one gets the XXX chain in a boundary 
magnetic field, which was discussed in the BA framework in [|13|] for < hi,h 2 < 2. From 
®~® one has in this limit 

h = rh- (35) 

There is only one critical field h cr = 2, which is the limit of h cr . Passing from summation 
to integration in eq. (PTf), we obtain for < hi,—h 2 < 00: 

1 1 7r r , e~ { ^- 1)x - e { ^- 1)x + e~ x 



E SU rf = -(hi + h 2 ) - - + - - / dx 



2 2 2 ./n cosh a; 
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where the second line was obtained from the first one after a simple manipulation. We can 



compare this result with that of |T3[ and find the agreement. For hi = hi = one has from 
([§): £ sw/ = (7T-l)/2-m2. 

Another aspect is the structure of the ground state in the regime h\^ > 0. Assuming 

(2) 

that, for example, for h\p > h\ r the ground state contains both left and right boundary 
1-strings to minimize the energy, we end up after a short calculation with a half-integer spin 
of the vacuum, which signals that such a state cannot, in fact, be the vacuum. Hence, the 
ground state must have a more intricate structure. Appealing to the Ising limit 7 — > 00, one 
sees that for > A the ground state must have both boundary spins directed opposite to 
the magnetic field, and therefore contain a kink in the bulk (recall that L is even). On the 
other hand, for h 12 < A the lowest energy configuration is such that the boundary spins are 
antiparallel, which means that the physical vacuum contains what was called a boundary 
excitation at one of the ends. This suggests that for finite A the correct ground state wave- 
function of the Hamiltonian ([I]) should contain a bulk hole with the minimal possible energy 
(i.e. the kink with zero rapidity 9 = 0) and both boundary 1-strings when hi 2 > h cr ■ Such 
a state has spin zero. Since the energy difference between this ground state and a state with 
a moving (instead of statinary) kink vanishes when the rapidity of the kink goes to zero, the 
spectrum acquires a gapless branch. Similarly, when hcr < ^1,2 < hcr , for the ground 
state to have the integer charge, it should also contain a kink in the bulk. When h\ t 2 < hcr 1 
the physical vacuum contains only one of the two boundary 1-strings and no stationary kink 
in the bulk (when hi = h 2 there are two possibilities to have either left or right boundary 



3 In the Ising limit 7 — > 00 the energy of the kink is independent of rapidity, and therefore this 
branch degenerates to the vacuum. 
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1-string in the vacuum, corresponding to the obvious two- fold degeneracy of the Ising ground 
state in this case). Such a state has a smaller energy for hi 2 < hUr than the one with a hole 

(2) 

in the bulk and two boundary strings, while for /i 12 > ha- the state with the bulk hole is 
energetically preferable, since in this case e& > Eh (see Fig. 2 and 0). This situation is in 
some sense analogous to the case of the periodic XX Z chain with odd L, where the ground 
state contains a kink. Let us mention that according to the above discussion the surface 
energy in the case /ii )2 > hu m is: 

E surf = -{hi + h 2 ) - g(A) + £ ft (0) - sinh 7 ( 1 + ^ ) . (37) 



n=l 

In the rational (7 — > 0) limit £^(0) vanishes and eq. (|37| ) becomes 



E surf = -(^1 + h) - - + - - / dx- ^— (38) 

2 2 2 y coshx 

This expression agrees with the one obtained in | 13fl . Note that the authors of |13| obtained 



eq. ( ^8l) under assumption that < < h cr , while our derivation shows that this result 
is valid for < hi 2 < 00. In the Ising limit eq. (|37|) gives the correct result E sur f = 
(3A -hi- h 2 )/2. 

In conclusion we would like to mention that within the BA technique it is possible to 
calculate also the boundary S-matrix for the scattering of kinks (represented by holes in 
the Dirac sea) in the ground state of the Hamiltonian (|l|) or in the excited boundary state. 
Such a calculation has been performed in |13[ for the boundary XXX chain and in ]I],|5| for 
the boundary sine-Gordon model. For the boundary XXZ chain these S-matrices have been 
obtained by Jimbo et al., |2j by the algebraic approach. 
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FIGURES 




FIG. 1. A schematic plot of the ground state surface energy as a function of the boundary 
magnetic field h = hi = — hi- The critical fields are also shown. 
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FIG. 2. Solid line: a schematic plot of the surface energy of the state with a boundary excitation 
at one of the ends, E sur f(h) + Sb(h), as a function of the boundary magnetic field h = h\ = —h.2- 
Dash-dotted line: the energy of the state containing a zero rapidity hole in the bulk E sur f(h)+Eh(0). 
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